Viral evolution under the pressure of an adaptive immune system - 
optimal mutation rates for viral escape 
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Based on a recent model of evolving viruses competing with an adapting immune system we 
study the conditions under which a viral quasispecies can maximize its growth rate. The range 
of mutation rates that allows viruses to thrive is limited from above due to genomic information 
deterioration, and from below by insufficient sequence diversity, which leads to a quick eradication of 
the virus by the immune system. The mutation rate that optimally balances these two requirements 
depends to first order on the ratio of the inverse of the virus' growth rate and the time the immune 
system needs to develop a specific answer to an antigen. We find that a virus is most viable if it 
generates exactly one mutation within the time it takes for the immune system to adapt to a new 
viral epitope. Experimental viral mutation rates, in particular for HIV (human immunodeficiency 
virus), seem to suggest that many viruses have achieved their optimal mutation rate. 



I. INTRODUCTION 



Since Eigen and Schuster introduced the concept of a 
quasispecies |^ , it has become a standard model to de- 
scribe molecular and viral evolution. If a simple, single- 
peaked fitness landscape is assumed, quasispecies the- 
ory predicts that error-prone replication leads to the for- 
mation of a central "master sequence" , surrounded by a 
cloud of mutant sequences. For viral evolution, this im- 
plies that any "wild-type" sequence is accompanied by a 
cloud of related mutants that, as a whole, represent a tar- 
get for the immune system. The quasispecies approach 
to molecular evolution has been the object of detailed in- 
vestigations, often supported by techniques of statistical 
physics I, |, ||, 0, |, I |1], |1], H g, @ revealing 
the characteristic features of such systems, including the 
occurrence of an error catastrophe. The latter charac- 
terizes a system in which a critical mutation rate exists 
beyond which the genomic information is irretrievably 
lost to mutations, i.e., beyond which selection ceases to 
operate [|ll], y[ |l6|, 0, |l|l (for an in-depth discus- 
sion of error catastrophes and related phenomena see also 
pO| ). The destabilizing effect of increased mutation rates 
has been observed for various viruses, including HIV pl| ] 
and Poliovirus p2[ . 

Recently, various extensions of the Eigen-Schustcr model 
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have been considered, in particular involving the shape 
of the fitness peaks and the landscapes' time-dependence. 
While the shape of the fitness function influences the ro- 
bustness of a species to mutations ^ , a behavior 
qualitatively different from the standard results can be 
observed for non- stationary fitness landscapes ^t\ . 
In rapidly changing environments, a second catastro- 
phe emerges besides the well-known error catastrophe, 
termed "adaptation catastrophe" . In a changing envi- 
ronment, sequence replication must occur with a non- 
vanishing error rate to enable the species to keep up with 
the environmental changes. (In static landscapes, a zero 
mutation rate is ultimately optimal because it maximizes 
the average global fitness of the population.) Indeed, 
a selective advantage for so called "mutator mutants" 
(or "general mutators" |2^) has been observed for Es- 
cherichia coli and Salmonella enterica under challenging 
living conditions @ 

For viruses in the environment of an adaptive immune 
system, the fitness landscapes for both the virus and 
the immune system are dynamically generated by a co- 
evolutionary process. This dynamics can be studied 
within the quasispecies' framework if the quasispecies 
character of both the viral population and the motifs of 
immune receptors is acknowledged. In an immune re- 
sponse, the presence of an antigenic epitope induces the 
proliferation of the corresponding immune receptor se- 
quence. This "master" sequence is associated with a 
cloud of closely related receptor sequences that emerge 
from somatic hypcrmutation of B-ccUs in the germinal 
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centers [ p2| . Competition between a viral population and 
an adaptive immune system takes place via an asymmet- 
ric coupling: while the immune quasispecies is strongly 
attracted by the virus, the viral quasispecies is driven 
away from its current master sequence by the immune 
system. This predator-prey-dynamics results in a migra- 
tion through sequence space as observed in many infec- 
tious diseases, such as HIV [|3[ 

The co-evolutionary dynamics within an infected host 
was recently formalized within a model relying only on 
a few dynamical rules recapitulated in the following 
section. Here, we focus on the implications of an optimal 
immune response within this framework, and consider 
the conditions that correspond to maximal viral fitness. 
Finally, we compare known viral mutation rates to those 
expected if a viral population has achieved an optimal 
mutation rate. 



II. VIRUS-IMMUNE SYSTEM CO-EVOLUTION 

Let us assume that the viral and the immunological qua- 
sispecies alike experience a single-peaked fitness function 
(Figure |l]), albeit one that can change in time. Let us as- 
sume further that at any particular time, the (viral) mas- 
ter sequence of length n grows at a rate tT„ (much larger 
than the "off-peak" or background- fitness rjy), and sim- 
ilarly for the immune system: ats ^ rjis. Such a simple 
immunological fitness function results from a reduction of 
the viral impact to induce proliferation of immune cells to 
its master sequence. Analogously, only the dominant im- 
mune sequence imposes a decay rate 5 on its complemen- 
tary viral sequence. Both species replicate imperfectly, 
with copy fidelities < I and qis < 1 (denoting the 
probability for correct duplication of a monomer drawn 
from an alphabet of size A). 




space: 

1. Once the immune system imposes a decay rate 
5 > on the viral master sequence (centered at the 
viral fitness peak), the narrow niche of the virus is 
assumed to move to an arbitrary sequence of the 
first error class. 

2. The viral quasispecies adapts to this new fitness 
peak on a time scale Ty . 

3. The fitness peak of the immune quasispecies is ad- 
justed, and moves to the new maximum of the viral 
distribution. 

4. The immune system adapts to the new fitness peak 
on the time scale Tis- 

As discussed previously the dynamically generated 
time scale t = + Tis can be approximated by the two 
expressions 
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The relative growth of the (moving) viral and immuno- 
logical master sequences in comparison to the environ- 
mental (background) sequences' growth can be deter- 
mined as Ip], pOf: 



{X - l){ai - T]i)qi 
leading to the conditions 



, ie{v,is}, 
(3) 

(4) 



for viability of the viral and immunological species, re- 
spectively. The regimes of (co-)existence of the two qua- 
sispecies can be determined by analyzing and Ki^. In 
particular, the viral quasispecies is subject to both a clas- 
sical error catastrophe at high mutation rates, and an 
adaptation catastrophe for small mutation rates. In con- 
trast, the immune system (as the driving force) is not 
subject to a limiting migration velocity, and accordingly 
only displays the classical error catastrophe . 



FIG. 1: Co-evolution of viral and immune quasispecies. 



III. OPTIMAL VIRAL MUTATION RATE 



The virus-immune system interaction is implemented by 
the following dynamic rules that are cyclicly iterated, 
leading to the quasispecies' migration through sequence 



Having derived the relations quantifying viral as well 
as immunological viability, we can now deduce optimal 
strategies for both the virus and the immune system. 



3 



The immune system attempts to minimize viral growth 
(^q^ = 0) which impHes the relation 
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between the optimal immune receptor size n^s and the 
per-site mutation probability His- This prediction and 
how it fares against the background of experimental data 
has been discussed in . Below, we extend this approach 
to derive the conditions for optimal viral escape from an 
immune response. 

Let us first approximate k„ in Eq. (^) by 



(6) 



using ay rjy, w 1. Optimizing viral viability con- 
ditions is akin to maximizing the viral species' relative 
growth rate Ky such that 
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Inserting r 
dition 



Tis into (p|) leads the equivalent con- 
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We can simplify this expression in the following man- 
ner. Writing \^ in terms of the mutation probability 
/ii, = 1 ~Qy rather than the copy-fidelity qy allows us to 
expand in terms of fiy (while leaving the term in In fiy 
untouched). Assuming furthermore that 5 ':S> cfy rjy, 
and n ^ 1, we find 



dhiy 



= <^ 6^ + n5cTy{hi ^y — 5Tis)^Jiy ~ Q ■ (9) 



We now proceed to determining the root of this expres- 
sion. While this can be done numerically (see below) , we 
first attempt to obtain an analytical approximation that 
permits an intuitive interpretation. For this purpose, it 
is allowable to assume In/i^, w const, as ln/i„ is a slowly 
varying function of fiy. The optimal per-site mutation 
probability /u* then follows as 
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Figure |^ shows a comparison between the optimal muta- 
tion rate /x* as given by the approximation (|K), and the 
exact solution /^* obtained numerically from (g). Despite 
the many approximations that have entered the deriva- 
tion of ( [Tc| ) , the analytic approximation is in good agree- 
ment with the numerical results. Improvements to the 
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FIG. 2: Optimal per-site mutation rate fi*, comparison be- 
tween the analytic approximation as given by equation jlC ) 
(dashed lines) and the numerical solution to equation (| ) 
(solid lines). Parameters are cr„ = 10, r^t, = 1, ais — 10, 
rjia — 1, Qis = 0.99, n = 100, 5 — 200, A = 4, unless specified 
otherwise in the plot. 



analytic approximation are possible if we neglect fewer 
of the higher order terms pl| . 

Let us now rewrite (^|) in terms of the optimal genomic 
mutation rate 
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This form suggests the following intuitive interpretation. 
The immune system adapts to a new virus strain within 
a time-span r^s , while the virus replicates in a time-span 
1/ay. The ratio between these two time scales measures 
the duration of one generation of the virus in units of the 
response time of the immune system. Hence, Equation 
([Tl) implies that the virus can optimally evade the im- 
mune system if the virus suffers on average one mutation 
per genome within the time the immune system needs to 
adapt to a new strain (Fig. ||). 

This condition guarantees that a maximal number of viri- 
ons have mutated away from the epitope to populate its 
first error class, precisely at that point in time when the 
immune system has adapted to attack the new viral qua- 
sispecies. 

If a viral quasispecies optimizes its mutation rate accord- 
ing to Eq. ( pT| ) , we expect to see this reflected in a relation 
between the mutation rate and genome size, such that 
their product is constant (given a particular generation 
time \/ijy). Optimization of genomic mutation rate can 
take place via an optimization of sequence length, given 
any particular per-site mutation rate. Table ^ shows that 
the genomic mutation rate only slightly varies within 
the class of RNA viruses, which presumably have a sim- 
ilar generation time. This is well in agreement with the 
prediction (pi]). 

Given the adaptation time of the immune system r^g 
and the generation time 1/cr^, we can test the predic- 
tion Eq. ( |ll| ) more specifically. The adaptation time Tig 
is the time necessary for the immune system to develop 
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TABLE II: Comparison of the genomic mutation rate of 
HIV-1 with the theoretical estimate (avTia)~^ from formula 
(pi]). Data are reproduced from [ ^o] , |4o[ 



Table || shows that the optimal genomic mutation rate 
as predicted by formula (O) is well within the range of 
the experimentally determined rate. This suggests that 
HIV-1 has adapted its mutation rate to optimally escape 
the immune system as well as the error catastrophe. 
IV. SUMMARY 



Tis ► 

FIG. 3: Regrowth from a single virus particle to a popu- 
lation size of eight, within the time-span r^s (dots indicate 
mutations). The virus can best evade the immune system if 
almost every virion in the population ai t = na differs from 
the initial virion by exactly one mutation. 



Organism ji /^i, = n/ii, 


Lytic RNA-bascd viruses [ 36 


1 


Poliovirus 


7.4 • 10^ 


1.1 ■ 10"'' 


0.81 


Influenza A Virus 


13.6 • 10-' 


> 7.3 ■ 10"^ 


0.99 


RNA-based Retroviruses |35 






Spleen Necrosis Virus 


7.8 • 10-'' 


2.0 ■ 10"^ 


0.16 


Molony Murine Leukemia Virus 


8.4 • 10^ 


> 3.5 ■ lO"** 


0.029 


Rous Sacroma Virus 


9.3 • 10^ 


4.6 ■ 10"^ 


0.43 


HIV-1 


9.2 • lO-'* 


2.4 ■ 10"^ 


0.22 



TABLE I: Genomic length n and spontaneous mutation rates 
per base pair and replication /i„ for RNA-based viruses that 
compete with advanced immune systems, as well as genomic 
mutation rate /^Jf = n^v Note that this product is an ap- 
proximation for — 1 — Ml))" for nfiy < 1. Data are 
reproduced from M, E^, 5j . 



The dynamics of co-evolution between virus and immune 
system can be studied within the framework of molec- 
ular evolution in time-dependent fitness landscapes, in 
which a constantly changing, polymorphic, viral popula- 
tion competes with an immune system adapting to keep 
track of the viral changes. Such an analysis reveals an 
optimal mutation rate for the immune system (so as to 
constrain the range of mutation rates within which the 
virus is stable) that appears to be compatible with avail- 
able data. The same formalism can be used to determine 
the optimal viral mutation rate, by maximizing the speed 
of adaptation while minimizing information loss due to 
mutations. It follows that the optimal viral mutation 
rate is reached if a sequence undergoes on average one 
mutation within the time it takes for the immune sys- 
tem to adapt to the viral genomic signature, thus barely 
staying ahead of the immune system. Such optimal muta- 
tion rates are compatible with experimentally determined 
ones, and suggest that the constancy of genomic muta- 
tion rates within viral classes (while sequence length and 
per-site mutation rates vary over many orders of magni- 
tude) can be explained by selection favoring viral strains 
at or near the optimal rate. 



a specific answer to an antigen. For most systems, this 
can be estimated to take between 7 to 14 days [Q. The 
generation times of viral species of course vary, but data 
from HIV-1 is available. 







r„ [d] 


(<T„r„)-i 




HIV-1 


0.4. ..3. 5 


7. ..14 


0.02. ..0.36 


0.22 
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